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Abstract We present world average values for D°-D° mixing parameters x and y, CP violation parameters 
\q/p\ and Arg(g/p), and strong phase differences S and S K7V7V . These values are calculated by the Heavy Flavor 
Averaging Group (HFAG) by performing a global fit to relevant experimental measurements. The results for 
x and y differ significantly from zero and are inconsistent with no mixing at the level of 6.7a. The results 
for \q/p\ and Arg(q/p) are consistent with no CP violation. The strong phase difference S is less than 45° at 
95% C.L. 
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1 Introduction 

Mixing in the D°-D° system has been searched 
for for more than two decades without success — un- 
til last year. Three experiments -Belle,' 1 ' Babar,' 2 ' 
and CDF [3 '- have now observed evidence for this 
phenomenon. These measurements can be combined 
with others to yield World Average (WA) values 
for the mixing parameters x = (m 1 — m 2 )/T and 
y = A(F]^ — r 2 )/(2r), where m 1 ,m 2 and r i; r 2 are 
the masses and decay widths for the mass eigenstates 
D 1 = p\D°)-q\D°) and D 2 = p\D°)+q\D°), and 
T = (T 1 +r 2 )/2. Here we use the phase convention 
CP\D°) = -\D°) and CP\D°) = -\D°). In the ab- 
sence of CP violation (CPV) , p = q = l/ \f2 and D 1 
is CP-even, D 2 is CP-odd. 

Such WA values have been calculated by the 
Heavy Flavor Averaging Group (HFAG)' 4 ' in two 
ways: (a) adding together three-dimensional log- 
likelihood functions obtained from various mea- 
surements for parameters x, y, and 8, where 8 
is the strong phase difference between amplitudes 
A(D° K+TT-) and A(D° -> K'TT+y, and (b) 
doing a global fit to measured observables for x, 
y, 8, an additional strong phase 8 Knw , and R D = 
\A(D ^K + TT-)/A(D ->K-7r+)\ 2 . For this fit, cor- 
relations among observables are accounted for by us- 
ing covariance matrices provided by the experimental 
collaborations. The first method has the advantage 



that non-Gaussian errors are accounted for, whereas 
the second method has the advantage that it is eas- 
ily expanded to allow for CPV. In this case three 
additional parameters are included in the fit: \q/p\, 
= Arg(?/p), and A D = (R+ - R~ D )/{R + D + R^), 
where the +(— ) superscript corresponds to D°(D a ) 
decays. When both methods are applied to the same 
set of observables, almost identical results are ob- 
tained. The observables used are from measurements 
of D° -> K+£~v, D° -> K+K-/-K+TT-, D° ->■ K+tt~ , 
D°^K + n-ir°, D^K+tt-jt+tt- , and D°^K° S tt+ti- 
decays, and from double-tagged branching fractions 
measured at the ij)(y!70) resonance. 

Mixing in heavy flavor systems such as that of 
B° and B° is governed by the short-distance box 
diagram. In the D° system, however, this diagram 
is doubly-Cabibbo-suppressed relative to amplitudes 
dominating the decay width, and it is also GIM- 
suppressed. Thus the short-distance mixing rate is 
tiny, and D°-D° mixing is expected to be dominated 
by long-distance processes. These are difficult to cal- 
culate reliably, and theoretical estimates for x and y 
range over two-three orders of magnitude.' 5, 6 ' 

With the exception of -0(3770) — > DD measure- 
ments, all methods identify the flavor of the D° 
or D° when produced by reconstructing the decay 
D*+^i^o 7r + Qr d*-^D°tt-; the charge of the ac- 
companying pion identifies the D flavor. For signal 
decays, M D ,—M D0 —M v+ = Q « 6 MeV, which is rela- 
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tively close to the threshold. Thus analyses typically 
require that the reconstructed Q be small to suppress 
backgrounds. For time-dependent measurements, the 
D° decay time is calculated via (l/p) x M D0 , where £ 
is the distance between the D* and D° decay vertices 
and p is the D° momentum. The D* vertex position 
is taken to be at the primary vertex' 3 ' (pp) or is cal- 
culated from the intersection of the D° momentum 
vector with the beamspot profile (e + e~). 

2 Input Observables 

The global fit determines central values and errors 
for eight underlying parameters using a \ 2 statistic 
constructed from 26 observables. The underlying pa- 
rameters are x, y, 5, R D ,A D , \q/p\, 4>, and S K7T7T . The 
parameters x and y govern mixing, and the param- 
eters A D , \q/p\, and (j) govern CPV. The parame- 
ter 5 Knn is the strong phase difference between the 
amplitude A(D° ^ K+ir-ir°) evaluated at M 



K+7 



if* (890) 



, and the amplitude A(D°^K 7r + 7r°) eval- 



uated at M„„ 

if 7 



■ M 



if* (890) ' 
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Fig. 1. WA value of R M from Ref. [4], as cal- 
culated from D°^K+£- v measurements. 1 

All input values are listed in Table 1. The ob- 
servable R M = (x 2 +y 2 )/2 measured in D° ^K + l~v 
decays' 7 ' is taken to be the WA value' 4 ' calculated by 
HFAG (see Fig. 1). The observables y CP and A r mea- 
sured in D°^K + K~ /n + Tr~ decays' 1, 81 are also taken 
to be their WA values' 4 ' (see Fig. 2). The observables 
from D°^KItt+tt- decays' 91 for no-CW are HFAG 
WA values,' 4 ' but for the CPV-allowed case only Belle 
values are available. The D° — > K + ir~ observables 
used are from Belle' 10 ' and Babar,' 2 ' as these measure- 
ments have much greater precision than previously 



published L>° -> K + Tr~ results. The D° -> K+tt't: 
and D" ^ K + tt~ tt + tt~ results are from Babar,' 111 and 
the V(3770)^L>£> results are from CLEOc.' 12 ' 

The relationships between the observables and the 
fitted parameters are listed in Table 2. For each 
set of correlated observables, we construct the dif- 
ference vector V, e.g., for D° — > A"g7r + 7r~ decays 
V = (Ax,Ay,A\q/p\,A(j)), where A represents the dif- 
ference between the measured value and the fitted pa- 
rameter value. The contribution of a set of measured 
observables to the \ 2 is calculated as V • (M _1 ) • V T , 
where M _1 is the inverse of the covariance matrix for 
the measurement. All covariance matrices used are 
listed in Table 1. 
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World average 
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H 
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3.420 ± 1.390 ± 0.740 < 



-1.200 ± 2.500 ± 1.400 1 
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Fig. 2. WA values of y CP (top) and A r (bot- 
tom) from Ref. [4] , as calculated from D° — » 
K + K~ /ir + Tr~ measurements.' 1 ' 81 
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3 Fit results 



The global fit uses MINUIT with the MIGRAD 
minimizer, and all errors are obtained from MINOS. 
Three separate fits are performed: (a) assuming CP 
conservation (A D and <p are fixed to zero, \q/p\ is 
fixed to one); (b) assuming no direct CPV (A D is 
fixed to zero); and (c) allowing full CPV (all param- 
eters floated). The results are listed in Table 3. For 
the CPV- allowed fit, individual contributions to the 
X 2 are listed in Table 4. The total x 2 is 23.5 for 
26 — 8 = 18 degrees of freedom; this corresponds to a 
confidence level of 0.17. 

Confidence contours in the two dimensions (x,y) 
or in (\q/p\,(f>) are obtained by letting, for any point 
in the two-dimensional plane, all other fitted param- 
eters take their preferred values. The resulting la-5a 
contours are shown in Fig. 3 for the CP-conserving 
case, and in Fig. 4 for the CPF-allowed case. The 
contours are determined from the increase of the x 2 
above the minimum value. One observes that the 
{x,y) contours for no-CP V and for CP^-allowed are 
almost identical. In both cases the x 2 a t the no- 
mixing point (x,y) = (0,0) is 49 units above the mini- 
mum value; this has a confidence level corresponding 
to 6.7(7. Thus, no mixing is excluded at this high 
level. In the (\q/p\,cj>) plot, the point (1,0) is on the 
boundary of the la contour; thus the data is consis- 
tent with no CPV. 



One-dimensional confidence curves for individual 
parameters are obtained by letting, for any value 
of the parameter, all other fitted parameters take 
their preferred values. The resulting functions Ax 2 = 
X 2 — Xmin (where xLin i s the minimum value) are shown 
in Fig. 5. The points where Ax 2 = 2.70 determine 
90% C.L. intervals for the parameters as shown in 
the figure. The points where Ax 2 = 3.84 determine 
95% C.L. intervals; these are listed in Table 3. 




Fig. 3. Two-dimensional contours for mixing 
parameters (x,y), for no CPV. 
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Fig. 4. Two-dimensional contours for parameters (x,t/) (left) and (\q/p\,(f>) (right), allowing for CPV. 



Chinese Physics C (HEP & NP) 



4 



Table 1. Input values used for the global fit, from Refs. [1, 2, 7-12]. 



Observable 


Value 


Comment 




VCP 
A T 


(1.132±0.266)% 
(0.123±0.248)% 


WA D°-+K+K-/ir 


+ 7r results [4] 


x (no CPV) 


(0.811 ±0.334)% 






y (no CPV) 


(0.309 ±0.281)% 


No CPV: 




\q/p\ (no direct CPV) 


0.95±0.22±°;J° 


WA D°^K° s tx+tt- 


results [4] 


4> (no direct CPV) 


(-0.035 ±0.19 ±0.09) rad 







\q/p\ 



(0.81±0.30±g;^)% 



(0.37±0.25±°;i°)% 
0.86±0.30±°;J° 
(-0.244±0.31±0.09) rad 



C*PV-allowed: 
Belle D°^K^tt+tt- 



1 

-0.007 
-0.255a 
0.216 



results. Correlation coefficients 
-0.007 -0.255a 0.216 
1 -0.019a -0.280 

-0.019a 1 -0.128a 

-0.280 -0.128a 1 



Note: a. = (\g/p\ ± 1) /2 is a variable transformation factor 



R M 


(0.0173 ±0.0387)% 


WA D°- 


^K+l~v results [4] 




x" 


(2.39±0.61±0.32)% 


Babar D°— >K^n n° result. Correlation coefficient 


= -0.34. 


y" 


(-0.14±0.60±0.40)% 


Note: x" = x cos 6 K7r7r + y sin 8 K7r7r , y" = y cos 5 Klv w - 


■xsm5 Knn . 


R M 


(0.019 ±0.0161)% 


Babar D ^K + tt~it + -k~ result. 








CLEOc results from "double-tagged" branching fractions 






measured in ip(3770)^DD decays. Correlation coefficients: 


R M 


(0.199±0.173±0.0)% 




1 -0.0644 0.0072 0.0607 




V 


(-5.207±5.571±2.737)% 




-0.0644 1 -0.3172 -0.8331 




R D 


(-2.395 ±1.739 ±0.938)% 




0.0072 -0.3172 1 0.3893 


> 


y/R D COS<5 


(8.878 ±3.369 ±1.579)% 




0.0607 -0.8331 0.3893 1 





Note: the only external input to these fit results are 
branching fractions. 



R 

x 



D 

12+ 



(0.303±0.0189)% 
(-0.024 ±0.052)% 
(0.98±0.78)% 




A 

x 
V 



D 
12- 



(-2.1±5.4)% 
(-0.020 ±0.050)% 
(0.96±0.75)% 



Babar D°^K + n results; correlation coefficients same as above. 



R 

x 



D 
12+ 



(0.364±0.018)% 
(0.032 ±0.037)% 
(-0.12±0.58)% 



Belle D°^K + n results. Correlation coefficients: 

1 0.655 -0.834 

0.655 1 -0.909 

-0.834 -0.909 1 





(2.3±4.7)% 




x' 2 - 


(0.006 ±0.034)% 


Belle D°^K + tt results; correlation coefficients same as above. 


y'~ 


(0.20±0.54)% 
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Table 2. Left: decay modes used to determine fitted parameters x, y, 8, 5 Kn7T , R D , A D , \q/p\, and <j>. Middle: 
the observables measured for each decay mode. Right: the relationships between the observables measured 
and the fitted parameters. 



Decay Mode 


Observables 


Relationship 


D°^K+K- /tt+tt- 


VCP 
A r 


2 ycP = (\l/p\ + \p/l\)y cos^ - (\q/p\ 
2A r = (\q/p\-\p/q\)y cos <f> - (\q/p\- 


— \p/q\) xs'm<j) 
f \p/q\)xsm(p 


D°-*K%-k+-k- 


X 

y 
\q/p\ 
<t> 






R M 


R M = {x 2 W)n 


D°^K + TT-7T° 


x" 


x" = x cos S K7V7V +ysin8 Knn 




(Dalitz plot analysis) 


y" 


y" = ycos8 K7V7V -xsm8 K7V7v 




D°~>K+Tr-n+TT- 


R M 


R M = (x 2 +y 2 )/2 




R M 






"Double-taggcd" branching fractions 
measured in tp(3770)—*DD decays 


y 

R D 
R D cos 8 


R M = (x 2 +y 2 )/2 








R D = ( R D + R d)/ 2 








A D = (R+-R-)/(R++R-) 






R D> R D 

x' 2 +, x ,2 ~ 

y'+, y'~ 


x' = x cos 8 + y sin 8 

y' = y cos 8 — x sin 8 

A M = (| g /p| 4 -l)/(| 9 /p| 4 + l) 

x'± = [(l±A M )/(l T A M r/\x'cos4 

j/± = [(l±A M )/(l T A M )]VV cos ^ 


± y' sin <j>) 
=F x' sin<j>) 
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Fig. 5. The function A\ 2 = X 2 ~ Xmin f° r fitted parameters x, y, 8, S Knn , \q/p\, and <f>. The points where 
A\ 2 = 2.70 (denoted by the dashed horizontal line) determine a 90% C.L. interval. 
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Table 3. Results of the global fit for different assumptions concerning CPV. 



Parameter 


No CPV 


No direct CPV 


CPV-allowed 


CPV-allowed 95% C.L. 


x (%) 


0.98±°;|? 


Q7+ - 27 
u - a ' -0.29 


0.97 jig g 


0.39-1.48 


y(%) 


0.75 ±0.18 


0.78l° ;i| 


0.78^ 


0.41-1.13 


*(°) 


21.6t£« 


23.411^ 


21.9t»;« 


-6.3-44.6 


R D (%) 


0.335 ±0.009 


0.334 ±0.009 


0.335 ±0.009 


0.316-0.353 


A D (%) 






-2.2 ±2.5 


-7.10-2.67 


k/p| 




0.95 t°o\l 


0.86±S;i« 


0.59-1.23 


<M°) 




-a-7±S:l 


-9- 6 ±1:! 


-30.3-6.5 


<w (°) 


30-8 


o 9 r-+25.0 
-25.7 


32.4125-1 


-20.3-82.7 



Table 4. Individual contributions to the x 2 f° r the CPF-allowed fit. 



Observable 


x 2 Ex 2 


VCP 
A r 


2.06 2.06 
0.10 2.16 


V K°TV + TV- 


0.20 2.36 
1.94 4.30 
0.00 4.30 
0.46 4.76 


R M ( K+e -~ v ) 


0.06 4.83 


y K+TT-Tt° 


1.24 6.06 
1.62 7.69 


R M /y/R D /y r R~^cos5 (CLEOc) 


5.59 13.28 


Rt/x' 2+ /y' + (Babar) 
R' D /x' 2 -/y'- (Babar) 
R+/x' 2 +/y'+ (Belle) 
R'nlx 12 - ly<- (Belle) 


2.54 15.82 
1.75 17.57 
3.96 21.53 
1.43 22.95 


R M (K+n-n+n-) 


0.49 23.45 
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4 Conclusions 

From the global fit results listed in Table 3 and 
shown in Figs. 4 and 5, we conclude the following: 

• the experimental data consistently indicate that 
D" mesons undergo mixing. The no-mixing 
point x = y = is excluded at 6.7a. The param- 
eter x differs from zero by 3. Oct; the parameter 
y differs from zero by 4.1a. The effect is pre- 
sumably dominated by long-distance processes, 
which arc difficult to calculate. Thus unless 
\x\ 3> \y\ (see Ref. [5]), it may be difficult to 
identify new physics from mixing alone. 

• Since y CP is positive, the CP-even state is 
shorter-lived, as in the K°-K° system. How- 
ever, since x also appears to be positive, the 
CP-even state is heavier, unlike in the K°-K° 
system. 

• It appears difficult to accomodate a strong 
phase difference S larger than 45°. 

• There is no evidence yet for CPV in the D°-D" 
system. Observing CPV at the level of sensi- 
tivity of the current experiments would indicate 
new physics. 
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